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Abstract 

Recently, planthoppers outbreaks have intensified across Asia resulting in heavy rice yield losses. The problem has been 
widely reported as being induced by insecticides while other factors such as global warming that could be potential drivers 
have been neglected. Here, we speculate that global warming may increase outbreak risk of brown planthopper 
(Nilaparvata lugens Stal.). We present data that demonstrate the relationship between climate variables (air temperature and 
precipitation) and the abundance of brown planthopper (BPH) during 1998-2007. Data show that BPH has become 
significantly more abundant in April over the 10-year period, but our data do not indicate that this is due to a change in 
climate, as no significant time trends in temperature and precipitation could be demonstrated. The abundance of BPH 
varied considerably between months within a year which is attributed to seasonal factors, including the availability of 
suitable host plants. On the other hand, the variation within months is attributed to fluctuations in monthly temperature 
and precipitation among years. The effects of these weather variables on BPH abundance were analyzed statistically by a 
general linear model. The statistical model shows that the expected effect of increasing temperatures is ambiguous and 
interacts with the amount of rainfall. According to the model, months or areas characterized by a climate that is either cold 
and dry or hot and wet are likely to experience higher levels of BPH due to climate change, whereas other combinations of 
temperature and rainfall may reduce the abundance of BPH. The analysis indicates that global warming may have 
contributed to the recent outbreaks of BPH in some rice growing areas of Asia, and that the severity of such outbreaks is 
likely to increase if climate change exaggerates. Our study highlights the need to consider climate change when designing 
strategies to manage planthoppers outbreaks. 
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Introduction 

Outbreak frequency of Nilaparvata lugens (brown planthopper; 
BPH) has been increasing in Asian rice growing countries in recent 
years (2005-2012) [1]. These trends are widely linked to adverse 
effects on BPH natural enemies of the increased use of broad 
spectrum insecticides for control of a range of pests [2-11]. They 
are rarely linked to environmental changes such as global warming 
that can directly or indirecdy influence outbreak trends in rice 
growing regions [12]. Insecticide resistance is another reason for 
the increase. 

Rice is extremely important for human food supply, yet a 
comprehensive understanding of how climate change affects rice 
pests is still premature. As the global climate continues to change, 
there is a need to understand how planthoppers respond to low- 
frequency (interannual or longer period) climatic variability. Such 
low frequency climatic variability is of interest because of its 
significance for the understanding and prediction of protracted 
climatic anomalies [13]. Since ongoing climate change has a 
profound effect on insect pests and changes their pest status [14] 
and population dynamics [15], climate change probably contrib- 
utes to recent destructive outbreaks of brown planthopper (BPH) 
Nilaparvata lugens (Stal) (Homoptera: Delphacidae). BPH is one of 



the most serious pests of rice in both temperate and tropical 
regions of East and South Asia and has become especially 
problematic over the past few years [16]. They are active 
throughout the year in the tropics [17]. In Bangladesh, where 
rice is grown year-round, they are prevalent throughout the year 
and can produce 8-11 generations per a year [personal 
communication, Nur Ahmed] . 

BPH's life cycle comprises three distinct stages: egg, nymph and 
adult. Both nymphs and adults direcdy damage rice plants through 
sucking the cell sap from the base (stem) of the plants [18] and by 
transmitting viruses such as rice ragged stunt (RRSV) and rice 
grassy stunt (RGSV) [19-20] which cause severe losses. In 2005- 
2006, more than 485,000 hectares of rice in southern Vietnam 
were severely affected by viral diseases seemingly spread by BPH, 
resulting in the loss of 828,000 tons of rice valued at US$120 
million [21]. 

Recently, outbreaks of BPH have become ever more common 
in Asia where farmers experience devastating losses [1], and have 
threatened rice production in parts of Thailand, Philippines, 
Indonesia, India, Bangladesh, Malaysia and China in the years of 
2009-2010. Global warming may have been a contributing factor 
because steadily warmer autumns have occurred since the 1990s. 
Changed climatic conditions of particular interest are milder 
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winters and warmer summers as well as changing patterns of 
precipitation. The latter includes increased risk of both extreme 
precipitation and severe dry spells. These conditions are predicted 
by the IPCC [22] to become more frequent due to changes in the 
global climate. Overall, temperature is predicted to increase 1.5- 
4.5°C during the present century [23] and precipitation to rise 10 
to 15% because a warmer atmosphere holds more water [24—25]. 
These changes could profoundly affect the population dynamics 
and the status of insect pests of crops [24], [26]. At the same time, 
migrant species may also change status to become major rice pests. 
Thus, in Japan JV. lugens was dominant until 1984, Sogotula Jurcifera 
by 1993 and Cnaphalocrosis medinalis up to now [27]. 

Reported BPH outbreaks are mostiy attributed to the applica- 
tion of insecticides in 2008-2012 [2-11]. Other factors that may 
induce the problem and threaten rice production in major rice 
producing countries are gready neglected although temperature 
has already started to rise with concurrent changes in precipitation 
patterns. For an array of insects these phenomena may affect e.g. 
phenology and pest status. Analysis of trapping data from several 
decades indicates that climate change also causes a change in the 
phenology of insects [28]. It is suggested that the use of the lower 
developmental threshold (To) and the thermal constant (A) that 
have been reported for insects would be useful in predicting the 
phenology of insect communities under global warming [29]. 
Therefore, a more detailed analysis of correlations between 
climatic parameters and phenology, abundance and pest status 
are needed. The hypothesis is that a higher frequency of 
destructive outbreaks of plant hoppers is related to climatic 
changes. 

In this paper, we use concurrent information on air tempera- 
tures, precipitation and abundance of BPH, recorded in 
Bangladesh over a 10-years period, to analyze the complex 
relationship between the climate variables and the dynamics of 
BPH. The aim of the present study is to provide a better 
understanding of the role of global warming in eliciting destructive 
outbreaks of BPH during the recent years as well as predicting its 
future role. 

Materials and Methods 

Study Site 

This study was conducted at a field site belonging to Bangladesh 
Rice Research Institute (24°0'0"N, 90°25'48"E), which is about 
36 km north of Dhaka, the capital of Bangladesh. The site, which 
is located in a rice-growing area, has been used mainly for insect 
pest monitoring and forecasting since 1980. 

Data Collection 

The BPH population was monitored daily by a light trap (Fig. 1 ) 
which was placed close to rice fields. The trap used a 100-watt 
tungsten light bulb positioned at 2 m height. Lighting hours were 
set for 1 2 h (from evening to dawn) and a strip with dichlorvos (or 
Vapona) was put in each collecting cage to kill the insects when 
caught. A steel drum was placed below the cages to collect all 
insects including minute ones and to allow for their easy removal 
from the cages. The bulb and Vapona were changed when 
needed. The insects in the cages were removed every morning and 
preserved temporarily in paper bags dated daily before processing. 
Trapped insects were identified, counted and deposited in the 
Entomology Division of Bangladesh Rice Research Institute. 
There were no other lights near the trap. Monthly BPH catches 
from January to December from 1998 through 2007 were then 
calculated and used in this study (Appendix to manuscripts). 




Figure 1. The light trap used to collect brown planthopper 
(BPH) in this study. It was placed adjacent to a rice field belonging to 
Bangladesh Rice Research Institute (BRRI). 
doi:1 0.1 371/joumal.pone.0091 678.g001 

Climate Data Source 

Climate variables (daily maximum and minimum temperatures, 
and precipitation) covering the same interval and study area were 
provided by the Plant Physiology Division of Bangladesh Rice 
Research Institute (www.brri.gov.bd). 

Statistical Analysis 

Trends in climate and abundance of BPH from 1998 
through 2007. We analyzed whether the climate data showed 
any long-term trends with respect to monthly temperatures and 
precipitation. For this purpose, we applied three climatic variables 
for every month through the 10 years period: MinTemp and 
MaxTemp which represent the average of the daily minimum and 
maximum temperatures during a month, and Rain which is the 
average amount of rain falling per day during the month. The 
three dependent variables were analyzed by a general linear model 
(PROC GLM in SAS [30]) with Month as a classification variable 
and Tear (number of years since 1998) as a quantitative variable. 
The model also included the interactions between Month and Tear 
to test whether trends differed among months. In the following, 
this model is termed the "null-model". Overall, the model 
included 12 month parameters, 1 year parameter, and 12 Month 
X Year interaction terms, to a total degree of freedom of 23. 

The monthly variation in BPH was analyzed in the same way as 
the climate variables, using the dependent variable jV calculated as 
the average number of adult BPH caught per day during a month. 
Prior to the analyses we checked the dependent variables for 
normality and variance homogeneity. Variables that did not meet 
these requirements were transformed with the appropriate 
transformation obtained by means of Taylor's power law [31]. 

Effect of weather variables on the abundance of 
BPH. The null-model predicting the expected number of BPH 
by means of the independent variables Month and Year and their 
interactions was expanded to include temperature and rainfall in 
order to investigate how much of the residual variance that could 
be accounted for by taking the actual weather in the month into 
account. The null-model expanded with the quantitative climate 
variables (called Temp and Rain) is termed the "full climatic 
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model". To avoid collinearity between the independent variables, 
minimum and maximum temperatures were not used in the same 
model because they tend to be positively correlated. 

The full climatic model included first and second order terms of 
Temp and Rain plus all interaction terms between these indepen- 
dent variables. Temp could be either MinTemp or MaxTemp. 
Variables were transformed if needed to meet the criteria of 
normality and variance homogeneity. Analyses were conducted by 
means of PROC GLMSELECT in SAS [30] using a backwards 
elimination procedure based on Akaike's corrected information 
criterion (AICC). The model with the lowest AICC was regarded 
as the best. 

Predicted effects of climate change on BPH 
population. We used the general linear model identified above 
to simulate what will happen to the BPH if climate change results 
in higher temperatures in combination with a changing pattern of 
rain. We applied the following scenarios: Daily average temper- 
atures were assumed to increase with 0°C, 1°C or 2°C and daily 
precipitation was either increased or decreased with 10%. The 
relative change in abundance of BPH for each month was 
calculated as R= 100(N-N 0 )/N 0 %, where M is the predicted 
catch of BPH for a given scenario and Mo is the corresponding 
value in the status quo scenario (i.e. the scenario with no climate 
change). 

Results and Discussion 

The period from December to February is the cold and dry 
season with low abundance of planthoppers (Fig. 2A-C). Insect 
pests often respond rapidly and dramatically to changes in climatic 
conditions affecting development (such as sudden precipitation 
and extreme temperatures), leading to large temporal variations in 
insect pest populations [32] . 

Trends in Climate and Abundance of BPH from 1998 
Through 2007 

Figure 2 indicates that standard deviations of rainfall and BPH 
tend to increase with the average. We therefore applied Taylor's 
power law [3 1] to identify the transformation of data that best 
stabilizes the variance. The appropriate transformations were 
obtained as y*=y 1 ' where p=\—b/2, b is the slope of the 
regression line fitted to the empirical values of log(variance) plotted 
against log(average) from each of the 12 months (Fig. 3). As p for 
daily rainfall is close to one third, the cubic root transformation 
was used to transform daily rainfall values into R — Rain 3 . For 
daily counts of PBH, b is close to 2, which means that this variable 
should be log transformed. As N in two cases was 0, we used the 
log(7V+l) transformation. There was no need of transforming 
temperature values. 

The model consisting of Year, Month and the interaction between 
them explained 97.2% of the variation in minimum temperatures, 
89.0% in maximum temperatures and 76.8% in the transformed 
values of daily precipitation. However, for all three variables, 
neither Year nor the interactions between Year and Month were 
significant (P>0.05). There was a considerable monthly variation 
in both temperatures and rainfall, but the analysis did not indicate 
that this pattern has changed over the 10-years period. 

With respect to a trend in BPH abundance, the null-model 
explained 84.8% of the variation in j* = log(jV+l) (^23,119 = 23.21; 
P<0.0001). There was no overall trend over years (P — 0.2468), 
whereas both Month and the interactions between Month and Year 
were significant {Month: P<0.0001; Year*Month: P = 0.0020). Only 
April showed a significant trend in BPH over the 10 years with a 
slope of 0. 1540 per year (P<0.0057). Month alone explained 79.4% 
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Figure 2. Monthly averages (±SD) of (A) Minimum and 
maximum temperatures; (B) daily rainfall; and (C) daily catches 
of Brown Planthoppers (BPH). Data from Bangladesh Rice Research 
Institute (BRRI). 

doi:10.1 371/journal.pone.0091 678.g002 

of the total variation in log (JV+1). As we could not demonstrate 
that the climate had changed, the most likely explanation for the 
increased abundance of BPH in April is that this month coincides 
with the growing and harvesting period of Boro rice, and that 
areas planted with Boro rice have expanded during the recent 
years [33-34]. 

Effect of Weather Variables on the Abundance of BPH 

When the null-model used above to predict the abundance of 
BPH was extended by including climatic variables, MinTemp had a 
higher predictive power than MaxTemp. We therefore used 
MinTemp to represent temperature in the following analyses. The 
full climatic model explained 86.3% of the variation in the 
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Figure 3. Log (variance) plotted against log(average) for (a) daily rainfall and (b) daily catches of BPH adults. Each dot represents a 
month. The straight line for daily rainfall is described by y= 1 .3068x+0.1252 (fi 2 = 0.9587) and for daily catches of BPH by y= 2.051 2x-0.2258 
(fi 2 = 0.9793). 

doi:1 0.1 371 /journal.pone.0091 678.g003 



observed values of log(jV+l) (i% 1>88 = 1 7.9; P<0.0001). The AICC 
of the null-model (i.e. the model without climate variables) was — 
94.80. The climate-driven model with the lowest value of AICC 
(-107.12) contained the terms Month, T, f\ R, R 2 , and f z R 2 , 
where Tis the monthly minimum temperature (i.e. T=MinTemp) 
and. R = Rain 1 ''. This reduced version of the full model (called the 
best climatic model) explained 83.6% of the variation in log(JV+l) 
(^16,103 = 32.8; P<0.0001). The parameters associated with the 
model's variables are shown in Table 1. 



Though the best climatic model explained slightly less of the 
variation than the null-model did (8.3.6% vs 84.8%), the former has 
7 parameters fewer than the null-model. However, the main 
advantage of the climate-driven model is that it provides a tool for 
analyzing how temperature and rain influence the dynamics of 
BPH. Thus, the model's parameters indicate that increasing 
temperatures will have a positive influence on BPH at low 
temperatures but a negative influence at high temperatures. Low 
levels of precipitation (e.g. <5 mm) will benefit BPH as long as the 
minimum temperature is below ca 23°C. Above this threshold, 
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Table 1. Parameter values of the model that best describes the observed values of log(A/+1), where N is the average number of 
BPH caught per day during a month. 





Parameter 


Estimate 


Standard Error 


t Value 


P 


January 


-2.3351 


1 .5228 


-1.53 


0.1282 


February 


-2.3138 


1 .7269 


-1.34 


0.1832 


March 


-1.6957 


1.8130 


-0.94 


0.3518 


April 


-0.2333 


1.8110 


-0.13 


0.8978 


May 


0.1209 


1.7931 


0.07 


0.9464 


June 


-1.0147 


1.7756 


-0.57 


0.5689 


July 


- 1 .2602 


1.7753 


-0.71 


0.4794 


August 


-0.6212 


1.7711 


-0.35 


0.7265 


September 


-0.5186 


1 .7802 


-0.29 


0.7714 


October 


0.2358 


1 .8074 


0.13 


0.8964 


November 


0.0983 


1.8015 


0.05 


0.9566 


December 


- 1 .4999 


1.6344 


-0.92 


0.3609 


T 


0.3528 


0.1767 


2.00 


0.0485 


7= 


-0.0120 


0.0044 


-2.69 


0.0083 


R 


0.7169 


0.2094 


3.42 


0.0009 


R 2 


-0.9471 


0.2213 


-4.28 


<.0001 


7* A* 


0.00119 


0.00029 


4.09 


<.0001 



T is the temperature and R the cubic root of the daily rainfall in month f. The predicted value of log(N+1) for month / is obtained as 

log(Af+l) = Af,+0.3528r-0.0120r 2 + 0.7169i?-0.9471^ 2 + 0.00119r 2 « 2 where M, is the parameter associated with month /'. P is the probability that the true 

parameter value is equal to 0. 

doi:1 0.1 371 /journal.pone.0091 678.t001 



more rain will increase BPH. Figure 4 shows the predicted 
combined effects of temperature and rainfall when added to the 
monthly value of log(jV+l). 

Predicted Effects of Climate Change on BPH Population 

Figure 5 shows the predicted monthly abundances of BPH 
under the various scenarios while Figure 6 shows the predicted 
relative changes in BPH abundance. The model predicts that BPH 



will become more abundant in January and during June-July, in 
particular if higher temperatures are associated with more 
precipitation. As these months are usually characterized by low 
abundances of BPH, the economic consequences are likely to be 
negligible. On the other hand, climate change is expected to 
decrease abundance of BPH during the spring and autumn 
seasons. These changes in seasonal patterns seem to benefit Boro 
rice crops (April-May), whereas rice grown in the Aus July- 
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Figure 4. The effects of minimum temperature and rainfall on the predicted abundance of BPH (log(/V<-1 )) (see Table 1 for further 
explanation). 

doi:1 0.1 371 /journal.pone.0091 678.g004 
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August) and Aman (November-December) seasons will become 
more exposed to BPH. Overall, the daily mean number of BPH on 
a yearly basis is predicted to decline with an increase in 
temperature (Fig. 7). The only scenario that may increase BPH 
is the one where precipitation declines without a concurrent 
increase in temperature. 

Even though our model indicates that BPH should become less 
abundant in Bangladesh if climate change causes an increase in 
the mean temperatures, it does not imply that this will be the case 
everywhere. Areas characterized by climatic conditions similar to 
those prevalent in Bangladesh during January (cool and dry) or 



July (hot and wet) seems most vulnerable to a changing climate. 
Using the model in combination with GIS-based climate 
information could assist in pointing out regions where future rice 
production is jeopardized. 

The results confirm that the combined effect of temperature and 
rainfall is complex [35-36]. Global warming may increase 
evaporation and precipitation with 10 to 15% because a warmer 
atmosphere holds more moisture [24-25]. In warm areas this is 
expected to increase populations of BPH. In some areas, however, 
global warming is predicted to result in a drier climate. Thus, Dai 
[37] suggests severe and widespread droughts in the next 30-90 
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Figure 5. The predicted daily catches of adult BPH (±SD) if minimum temperature increases with either 0°C, 1 °C or 2°C, and daily 
precipitation either decreases or increases with 10%. 

doi:1 0.1 371 /journal.pone.0091 678.g005 
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Figure 6. The predicted relative changes in daily catches of BPH under the various scenarios. 

doi:1 0.1 371 /journal.pone.0091 678.g006 



years over many land areas. According to the model, this may also 
exaggerate problems with BPH, but only in areas where minimum 
temperatures stay below ca 23°C. 

Climate change, and particularly global warming, may have a 
dramatic impact on pest insect species [38] . Many authors predict 
that climate change will have a number of effects on insects: 
sweeping shifts in herbivory rates; altered distribution and 
outbreak frequency of key insect pests; unpredictably altered 
relationships with natural enemies; and a general decrease in 
biodiversity [39-41]. The findings of the present study indicate 
that the reported BPH outbreak in rice fields occurring during the 
warm seasons might be attributed to global warming [42]. 



However, the model also shows that temperatures exceeding 
23°C in a dry climate may have an inhibiting effect on population 
density, implying that the consequences of global warming on 
planthoppers are ambiguous. Thus, their long-term abundance 
may increase during the cold season and decrease during the hot 
season, depending on how the ambient temperatures change 
relative to the threshold value of ca 23°C. Planthopper outbreaks 
were observed in the boro rice growing season i.e. the dry season 
[43] and significant BPH outbreaks were recorded during the dry 
season in the early 1990s in the Vientiane Plain [44] along with 
higher temperatures. The population subsequently declined in 
January due to low temperatures and peaked again in March and 
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Figure 7. The predicted mean daily catches of BPH averaged over a year if daily minimum temperature changes with either 0°C, 
+1 °C or +2°C, and if daily rainfall changes with either 0%, -10% or +10%. 

doi:1 0.1 371 /journal.pone.0091 678.g007 



April 2008 during the dry season [45] when temperatures were 
high. These weather conditions seem to promote BPH outbreaks, 
which is in accordance with the predictions based on our model, 
because minimum temperatures in the area during the dry season 
rarely exceed 23°C [12]. 

Wide-scale outbreaks of BPH in tropical Asian rice in the 1970s 
and 1980s [46-47] were attributed to the insecticidal destruction 
of natural enemies [48]. One of the main causes of the recent 
series of outbreaks in Thailand and Indonesia was also insecticide 
misuse [49]. Although insecticide applications to rice do not 
always trigger BPH outbreaks [46], they often disrupt the actions 
of BPH's natural enemies either by their direct killing action or by 
disrupting food chains [50] . Spraying of insecticide in rice field for 
controlling other pests destroys essential ecosystem services that 
regulate invading planthopper adults, thus increasing the risk of 
hopperburn on the farm [1 1] . Sublethal insecticide applications 
could theoretically increase BPH's capacity for migration because 
the planthoppers would acquire more fat and sugar, which provide 
fuel for flight, when feeding on insecticide-treated rice plants than 
when feeding on untreated plants [49]. The applications of 
sublethal doses of certain insecticides could enhance both BPH's 
reproductive and migratory capacity and theoretically increase the 
threat of BPH outbreaks even if the insecticides do not harm 
natural enemies [51]. 

It is tempting to conclude that human-controlled inputs such as 
insecticides or insecticides in combination with nitrogen fertilizer, 
for example, have been totally responsible for the synchronous 
BPH outbreaks in tropical rice observed across many areas of Asia 
shortly after the beginning of the green revolution and again more 
recently [49], but there is a widespread concern that outbreaks of 
herbivorous insects will increase in frequency and severity as a 
consequence of global warming. These predictions, however, are 
untested and still very uncertain. 

In an environment conducive to more rapid BPH increases, 
natural enemies and BPH-resistant cultivars would be expected to 
have less effect in regulating the density of BPH, and pesticides 
harming natural enemies would therefore be expected to have an 
above-average negative impact [49]. Kiritani [27] also reported 



that every 1°C rise from around X— 4°C would result in a 
decrease in winter mortality of about 16.5% based on a regression 
model linking winter mortality (T) to the mean temperature in 
January IX). 

jV. lugens is a long-range migrant species that immigrates into 
Japanese paddy fields from mainland China in early summer [52]. 
Our study does not consider immigration of BPH into the study 
area because this phenomenon has not yet been observed in 
Bangladesh. Moreover, we accentuate that our model is a simple 
measure of the importance of climate compared to all other factors 
(such as pesticides, natural enemies etc). It does not focus on the 
overall causes for pest outbreaks, nor does it independendy 
attribute pest trends to the many technological, climatic and biotic 
factors that may affect such trends. However, it discloses a 
pragmatic approach to assessing the combined importance of 
climate variables with special reference to evaluating the possible 
impact of climate change on pest outbreaks. 

Long time series like ours provide a good basis for linking pest 
outbreaks to climate change and to separate the influence of other 
factors, such as pesticides, from that of climate. Our study 
indicates that climate change may have contributed to the recent 
destructive outbreaks of planthoppers in rice crops [53]. However, 
more research is needed to predict more accurately the long-term 
consequences of a changing climate, and especially how such 
changes will affect pest dynamics at both regional and local scales. 

Conclusion 

Data collected from 1998 through 2007 did not indicate that 
temperature and precipitation patterns have changed significandy 
during this 10-years period. Therefore, we attribute the significant 
increase in BPH during April to changes in agricultural practice 
rather than climate change. However, the statistical models 
demonstrate that temperature and rainfall are driving variables 
for the dynamics of BPH, so any changes in these factors due to 
climate change are likely to affect the outbreak patterns of BPH. 
We tested this by simulating different climatic scenarios and found 
a quite complex picture depending on the severity and direction of 
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the expected climatic changes. Thus, our study emphasizes the 
importance of taking future climate into consideration when 
deciding which, where and when the different rice cultivars should 
be planted and harvested so as to minimize crop losses due to 
pests. 

Given the challenge of linking pest impacts and directional 
climate change for well-studied agricultural, maricultural (cultiva- 
tion of marine organisms in their natural habitats, usually for 
commercial purposes), and human diseases, it is not yet possible to 
predict the consequences for biodiversity. Very few empirical 
studies directly explore the relationship between climate and rice 
pest problems. Even fewer explore interactions between temper- 
ature and various components of an insect's life cycle [54-56]. We 
therefore identify three priorities for research to improve our 
ability to predict impacts of climate change on rice pests: 

1. Collect baseline data on insects of both major and minor importance: 
Baseline data are critical to predict changes in a warmer 
climate, but such data are rarely collected for rice ecosystems. 
Monitoring programs for the prevalence and severity of pests 
and their population- and community-level impacts must be 
implemented for a wider range of natural systems. 

2. Separate the effects of multiple climate variables on insects: To predict 
accurately future responses to climate change, we must 
quantify the direct and synergistic effects of multiple climate 
variables, such as temperature and precipitation, on insects. 
Laboratory and/or field experiments are crucial for studying 
the specific effects of these variables [57]. 

3. Forecast epidemics: Forecasting models using climate variables can 
effectively predict outbreaks for some crops and pests. Crop 
disease programs have long been in effect as, e.g., rice blast 

References 

1. Heong KL (2010) Assessing BPH outbreak risks of commonly used insecticides. 
Available: http:/ /ricehoppcrs. net/20 10/09/farmers-in-central-thailand-remain- 
trappcd-by-thc-bph-problcm/. Accessed 30 Sept 2010. 

2. Chien HV, Cuong L (2009) Farms that apply insecticides for leaf folder control 
are 10 times more at risk to hopperburn. Available: http://ricehoppers.net/ 
2009/09/ farms-that - apply- insecticide s-for-leaf-f blder-control-are- 1 0-timcs- 
more-at-risk-to-hopperburn/. Accessed 13 Sept 2009. 

3. Heong KL (2009) Planthopper outbreaks in 2009. Available: http:/ / ricehoppcrs. 
net/2009/09/. Accessed 25 Sept 2009. 

4. Islam Z, Haque SS (2009) Rice planthopper outbreaks in Bangladesh. Available: 
http:/ / ricehoppers.net/2009/08/ rice-planthopper-outbrcaks-in-bangladesh/ 
planthopper-outbreaks-in-2009. Accessed 20 Aug 2009. 

5. Moni (2011) BPH continues to threaten Thai rice farmers — heavy losses 
expected. Available: http:/ / ricchoppers.net/201 1 /04/bph-continucs-to- 
threaten-thai-rice-farmers heavy- losses-expected/. Accessed 20 April 2011. 

6. Moni (2012) Planthopper problems intensify in Thailand's rice bowl. Available: 
http://ricehoppers.net/2012/03/planthopper-problems-intensify-in-thailands- 
rice-bowl/. Accessed 21 March 2012. 

7. Luecha M (2010) Farmers' insecticide selections might have made their farms 
vulnerable to hopperburn in Chainat, Thailand. Available: http: //ricehoppcrs. 
net/2010/01/farmers%e2%80%99. Accessed 17 Jan 2010. 

8. Matsumura M, Sanada-Morimura S (2010) Recent status of insecticide 
resistance in Asian rice planthoppers. JARQ_ 44: 225-230. 

9. Soitong K, Sriratanasak W (2012) Will planthoppers continue to threaten 
Thailand's rice production in 2012? Rice Department, Bangkok, Thailand. 
Available: http://ricehoppers.net/20 1 2/02/will-planthopper-continue-to- 
threaten-thailands-ricc-production-in-2012/. Accessed 13 Feb 2012. 

10. Heong KL, Schoenly KG (1998) Impact of insecticides on herbivore -natural 
enemy communities in tropical rice ecosystems. In: Haskell PT, McEwcn P, 
editors. Eco toxicology: Pesticides and Beneficial Organisms. Chapman and Hall, 
London, UK. pp. 381-403. 

1 1. Teo G (201 1) Insecticide abuse in rice production causes planthopper outbreaks. 
Available: http://www.asianscicntist.com/topncws/irri-ban-insecticides-rice- 
production-duc-to-planthopper-outbrcaks-201 1/. Accessed 21 Dec 2011. 

12. BRRI (2012) Bangladesh Rice Research Institute annual report, 2012. 

13. Upmanu L, Michael M (1994) The Great Salt Lake: A barometer of low 
frequency climate variability. Water Resources Res 31: 2383-2635. 

14. Estay SA, Lima M, Labra FA (2009) Predicting insect pest status under climate 
change scenarios: combining experimental data and population dynamics 
modelling. J Appl Entomol 133: 491-499. 



(Pyricularia oryzae) models based on temperature and moisture 
forecasts when an epidemic will start and when to apply 
fungicide for optimal control [58]. Such forecasting programs 
are also in development for human diseases with climate 
sensitivity, such as Rift Valley fever, which is associated with 
warm El Nino events of high rainfall [58—59], and cholera, 
which is predictable from sea surface temperature associations 
with El Nino [59]. Similar forecasting models will be very 
useful for predicting epidemic outbreaks of insect pests in rice 
and we hope that the present paper will contribute to achieving 
this goal. 

Supporting Information 

Appendix SI Data used in this study. 

(XLS) 

Acknowledgments 

We greatly acknowledge the Bangladesh Rice Research Institute (www. 
brri.gov.bd) for supporting the staff to operate the data monitoring device 
including logistic support to collect field data, Entomology Division Staff 
for technical support, and Plant Physiology Division for supplying climate 
data. 

Author Contributions 

Conceived and designed the experiments: MFR NA MP A. Performed the 
experiments: MPA NA MAB. Analyzed the data: MPA GN DH. 
Contributed reagents/materials/analysis tools: MFR GN MPA. Wrote 
the paper: MPA GN DH. 



15. Karuppaiah V, Sujayanad GK (2012) Impact of climate change on population 
dynamics of insect pests. World J Agri Sci 8: 240—246. 

16. Hu G, Cheng XN, Qi GJ, Wang FY, Lu F, et al. (201 1) Rice planting systems, 
global warming and outbreaks of Mlapawata lugens (Sta.1). Bull Entom Res 101: 
187-99. 

17. KulshreshthaJP, Anjancyulu A, Padmanabhan SY (1974) The disastrous brown 
plant-hopper attack in Kerala. Indian Farming 24: 5-7. 

18. Ali MP, Alghamdi SS, Begum MA, Alam MZ, Huang D, et al. (2012) Screening 
of rice genotypes for resistance to the brown planthopper, Mlapawata lugens Stal. 
Cereal Res Commun 40: 502-508. 

19. Hibino H (1989) Insect-borne viruses in rice. In: Harris KF editor. Advances in 
disease vector research. Springcr-Verlag, New York, USA, pp. 209—241. 

20. Hibino H (1996) Biology and epidemiology of rice viruses. Annu Rev 
Phytopathol 34: 249-274. 

21. Du PV, Cabunagan RC, Cabauatan PQ, Choi HS, Choi IR, et al. (2007) 
Yellowing syndrome of rice: Etiology, current status, and future challenges. 
Omonrice 15: 94-101. 

22. IPCC (2007) Climate Change 2007. In: Solomon S, et al., editors. The Physical 
Science Basis. Contribution of Working Group I to the Fourth Assessment 
Report of the Intergovernmental Panel on Climate Change. Cambridge Univ. 
Press, Cambridge, UK and New York, NY, USA. 

23. Mcehl GA, Stocker TF, Collins W, Friedlingstein P, Gayc A, et al. (2007) Global 
climate projections. In: Solomon S, et al., editors. Climate Change 2007: The 
Physical Science Basis. Contribution of Working Group I to the Fourth 
Assessment Report of the Intergovernmental Panel on Climate Change. 
Cambridge University Press, Cambridge, pp. 747-845. 

24. Crosson P (1997) Impacts of climate change on agriculture, Climate Issues Brief 
No. 4. 

25. Cannon R (1998) The implications of predicted climate change for insect pests in 
the UK, with emphasis on non-indigenous species. Glob Change Biol 4: 785- 
796. 

26. Woiwod I (1997) Detecting the effects of climate changes on Lepidoptera. J Insect 
Conscrv 1: 149-158. 

27. Kiritani K (2006) Predicting impacts of global warming on population dynamics 
and distribution of arthropods in Japan. Popul Ecol 48: 5-12. 

28. Esbjcrg P (201 1) Phenology and pest status change of Agrotis segetum as related to 
the climate change in Northwestern Europe. IOBC/WPRS Bull 65: 43-49. 

29. Kiritani K (1988) Effects of climate change on the insect fauna (in Japanese). 
Mctcorol Res Rept 162: 137-141. 



PLOS ONE | www.plosone.org 



9 



March 2014 | Volume 9 | Issue 3 | e91678 



Warmer Climate Favours Outbreaks of Planthoppers 



30. SAS (2008) SAS Enterprise Guide 4.3. SAS Institute Inc. SAS Campus Drive, 
Gary, NC 27513. 

31. Taylor LR (1961) Aggregation, variance, and the mean. Nature 189: 732-735. 

32. Bale JS, Masters GJ, Hodkinson ID, Awmack C, Bczcmcr TM, ct al. (2002) 
Herbivory in global climate change research: direct effects of rising temperatures 
on insect herbivores. Glob Change Biol 8: 1—16. 

33. Singh UP (2002) Boro rice in Eastern India. Rice-wheat consortium regional 
technical coordination committee meeting. 10-14 February 2002. Rice-wheat 
consortium for the Indo-Gangctic Plains, New Delhi, India. Available http:// 
www.Scirp.org/journal/jep/. Accessed Nov 2009. 

34. Rahman MK, Parvin L (2009) Impact of irrigation on food security in 
Bangladesh for the past three decades. J Environ Protect 1 : 40-49. 

35. Escalada MM (2010) Thailand: Insect influx threatens rice production. 
Available: http://www.irinnews.org/report/90269. Accessed 25 Aug 2010. 

36. Petzoldt C, Seaman A (2008) Climate change effects on insects and pathogens. 
In: University of Vermont Extension, New York State College of Agriculture 
and Life Sciences, editors. Climate change and agriculture: promoting practical 
and profitable responses. New York State IPM Program, New York State 
Agricultural Extension Station, Geneva, New York. 

37. Dai A (2013) Increasing drought under global warming in observations and 
models. Nature Climate Change 3: 52-58. 

38. Logan JA, RegniereJ, Powell JA (2003) Assessing the impacts of global warming 
on forest pest dynamics. Front Ecol Environ 1: 130-137. 

39. Williams DW, Liebhold M (1995) Herbivorous insects and global change: 
potential changes in the spatial distribution of forest defoliators. J Biogeogr 22: 
665-71. 

40. Fleming RA, Volney WJA (1995) Effects of climate change on insect defoliator 
population processes in Canada's boreal forest: some plausible scenarios. Water 
Air Soil Poll 82: 445-54. 

41. Colcy PD (1998) Possible effects of climate change on plant /herbivore 
interactions in moist tropical forests. Climatic Change 39: 455-72. 

42. Schiller JM, Chaphcngxay MB, Linquist B, Rao SA (2006) Rice in Laos. Loss 
Banos (Philippine), IRRI. 435p. 

43. Win SS, Muhamad R, Ahma ZAM, Adam NA (201 1) Population fluctuations of 
brown planthoppcr. ,\/hipe//rat.a lugens Stal. and white backed planthopper, 
Sogatella j'urcifera Horvath on rice. J Entomol 8: 183-190. 

44. IRIN (2012) Humanitarian news and analysis. Thailand: Rice pests multiply 
post-floods. Available: http://www.irinnews.org/report/95058/Thailand-rice- 
pests-multiply-post-floods. Accessed 12 March 2012. 

45. Kontodimas DC, Eliopoulos PA, Stathas GJ, Economou LP (2004) Comparative 
temperature-dependent development of .Nephus includens (Kirsch) and Nephus 
hisignatus (Boheman) (Coleoptera: Coccinellidac) preying on Planococcus citri 



(Risso) (Homoptera: Pseudococcidae): Evaluation of a linear and various 
nonlinear models using specific criteria. Environ Entomol 33: 1-11. 

46. Heinrichs EA, Aquino GB, Chelliah S, Valencia SL, Reissig WH (1982) 
Resurgence of Mlapawata lugens (Stal) populations as influenced by method and 
timing of insecticide applications in lowland rice. Environ Entomol 11: 78-84. 

47. Shepard BM, Khan ZR, Pathak MD, Heinrichs EA (1991) Management of 
insect pests of rice in Asia. In: Pimentel D, editor. CRC Handbook of Pest 
Management in Agriculture. CRC Press, Boca Raton, FL, pp. 255—278. 

48. Matteson PC (2000) Insect pest management in tropical Asian irrigated rice. 
Annu Rev Entomol 45: 549-574. 

49. Bottrell DG, Schoenly KG (2012) Resurrecting the ghost of green revolutions 
past: The brown planthopper as a recurring threat to high-yielding rice 
production in tropical Asia. J Asia-Pacific Entomol 15: 122-140. 

50. Hcong KL, Schoenly KG (1998) Impact of insecticides on herbivore -natural 
enemy communities in tropical rice ecosystems. In: Haskell PT, McEwen P, 
editors. Ecotoxicology: Pesticides and Bcncfical Organisms. Chapman & Hall, 
London, UK, pp. 381-403. 

51. YinJL, Xu HW, Wu JC, Hu JH, Yang GQ, (2008) Cultivar and insecticide 
applications affect the physiological development of the brown planthopper, 
Mlapawata lugens (Stal) (Ilemiptera: Delphacidae). Environ Entomol 37: 206 212. 

52- Kiritani K (2001) Insects and cimate (in Japanese). Scizando, Tokyo. 

53. Hart AJ, Bale JS, Fcnlon JS (1997) Developmental threshold, day-degree 
requirements and voltinism of the aphid predator Episyrphus balteatus (Diptera: 
Syrphidae). Ann appl Biol 130: 427-437. 

54. JaramilloJ, Chabi-Olaye A, Kamonjo C, Jaramillo A, Vega FE, et al. (2009) 
Thermal tolerance of the coffee berry borer ILipulliein'mn^ hamper. Predictions of 
climate change impact on a tropical insect pest. Plos One 4(8): c6487. 
doi: 10. 137 1 /journal.pone.0006487. 

55. Harvell CD, Mitchell CE, WardJR, Altizer S, Dobson AP, et al. (2002) Climate 
warming and disease risks for terrestrial and marine biota. Science 296: 2 158— 
2162. 

56. Ishigua K, Hashimoto A (1991) Computer-based forecasting of rice blast 
epidemics in Japan. In: Teng PS, Polland LR, Argosino G, editors. Rice blast 
modeling and forecasting. International Rice Research Institute, Seoul, South 
Korea, pp. 39—53. 

57. Linthicum KJ, Anyamba A, Tucker GJ, Kelley PW, Myers MF, et al. (1999) 
Climate and satellite indicators to forecast rift Valley fever epidemics in Kenya. 
Science 285: 397. 

58. Anyamba A, Linthicum KJ, Tucker GJ (2001) Climate-disease connections: Rift 
valley fever in Kenya. Cad Saude Publica 17: 133-140. 

59. Pascual M, RodoX, Elmer SP, Colwell R, Bouna MJ (2000) Cholera dynamics 
and El Nino-Southern oscillation. Science 289: 1766. 



PLOS ONE | www.plosone.org 



10 



March 2014 | Volume 9 | Issue 3 | e91678 



